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The interplay between the finite volume and finite lattice spacing is investigated using lattice QCD 
simulations to compute the Landau gauge gluon propagator at zero temperature. Comparing 
several ensembles with different lattice spacings and physical volumes, we conclude that the 
dominant effects, in the infrared region, are associated with the use of a finite lattice spacing. The 
simulations show that decreasing the lattice spacing, while keeping the same physical volume, 
leads to an enhancement of the infrared gluon propagator. Moreover, we also present results for 
the Landau gauge gluon propagator at finite temperature. 
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1. The gluon propagator at zero temperature 

In [jl|], we have studied the lattice Landau gauge gluon propagator, using several ensembles 
with different lattice volumes and lattice spacings. Here we will report on the comparison of the 
gluon propagator computed from similar physical volumes but different lattice spacings. For details 
on the lattice setup see [[I]]. 

Figure |j] shows the gluon propagator, renormalized at 4 GeV, for physical volumes ranging 
from 3.3 fm to 8.1 fm. In all cases, for momenta above ~ 900 MeV the lattice data define a unique 
curve; this means that the renormalization procedure has been able to remove all dependence on 
the ultraviolet cutoff for the mid and high momentum regions. 

However, we observe in the infrared region sizeable effects due to the finite lattice spacing of 
the simulation. In particular, the plots show that large lattice spacing simulations underestimate 
the propagator in the infrared region. Moreover, we found that the corrections due to the finite 
lattice spacing are larger than the corrections due to the finite volume of the lattice simulation — 
for details see 

One can ask whether the analysis reported here depends on the renormalization scale. Usually, 
one chooses a renormalization point /J. in the high momentum region, in order to allow a contact 
with the perturbative behaviour. However, choosing a renormalization point in the infrared is also 
possible. Here we consider two alternative renormalization points, jU = 500 MeV and jU = 1 GeV 
and study how the propagators differ. Note that we only show results for the simulations performed 
at La 8 fm. 

Figures || and || show that for both renormalization points, the differences in the infrared 




Figure 1: Comparing the renormalized gluon propagator at jU = 4 GeV for various lattice spacings and 
similar physical volumes. 
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Figure 2: Gluon propagator renormalised at ji = 500 MeV (left) and at = 1 GeV (right). 



region are well beyond one standard deviation. This means that the dominant effect in the infrared 
propagator is indeed due to the finite lattice spacing, and not the physical box size of the simulation; 
this conclusion is independent of the renormalization scale. Therefore lattice artifacts due to a finite 
lattice spacing play an important role in the determination of the infrared behaviour of the gluon 
propagator in the Landau gauge. 





Figure 3: The gluon propagator renormalised at \x = 500 MeV (left) and at jj. = 1 GeV (right) at intermediate 
momentum range (top) and up to 5 GeV (bottom). 
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2. The gluon propagator at finite temperature 

Besides the Landau gauge gluon propagator at zero temperature, we are also engaged in the 
computation of the gluon propagator at finite temperature, in order to understand how temperature 
changes the gluon propagator [|2[ [|] . 

Finite temperature is introduced on the lattice through the reduction of the temporal lattice 
size: the simulations are done on lattices x L t , with L t <C L s . The temperature is defined by 
T = l/aL t . 

At finite temperature, the gluon propagator is described by two tensor structures, 

Df v (q) = 8 ah {Pj; v D T (q 4 ,q) + P^ v D L (q 4 ,q)) (2.1) 
where the transverse and longitudinal projectors are defined by 

P; v = (l-V)(l-M(s,v-^) , P^=(s MV -^)-P; v ; (2.2) 
the transverse Dj and longitudinal Di propagators are given by 

Dr(q) = 2v{N \_ l) ({A c !{q)A«(-q)) ~ ^(A^ (-<?)>) (2.3) 

= V{N \_ X) (i + f <^(<7M2(-<7)>) (2-4) 

The lattice setup is shown in table ||. The finite temperature simulations described in this 
section have been performed with the help of Chroma library [||] ; the FFT transforms have been 
done with the PFFT library [Q]. For the determination of the lattice spacing we fit the string tension 
data in [B], using the functional form used in [0], in order to have a function a(j8); the fit has a 
X 2 /dof ~ 0.03. 



Temp. (MeV) 


P 


L s 


U 


a [fm] 


1/a (GeV) 


121 


6.0000 


32,64 


16 


0.1016 


1.9426 


162 


6.0000 


32,64 


12 


0.1016 


1.9426 


243 


6.0000 


32,64 


8 


0.1016 


1.9426 


260 


6.0347 


68 


8 


0.09502 


2.0767 


265 


5.8876 


52 


6 


0.1243 


1.5881 


275 


6.0684 


72 


8 


0.08974 


2.1989 


285 


5.9266 


56 


6 


0.1154 


1.7103 


290 


6.1009 


76 


8 


0.08502 


2.3211 


305 


5.9640 


60 


6 


0.1077 


1.8324 


305 


6.1326 


80 


8 


0.08077 


2.4432 


324 


6.0000 


32,64 


6 


0.1016 


1.9426 


486 


6.0000 


32,64 


4 


0.1016 


1.9426 



Table 1: Lattice setup used for the computation of the gluon propagator at finite temperature. 
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■— • 32 3 xl6 0=6.0 (T = 121 MeV) 
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— ■ 32 3 xl4 0=6.0 (T = 139 MeV) 




*— « 32 3 xl2 0=6.0 (T = 162 MeV) 




— * 32"'xl0 0=6.0 (T = 194 MeV) 


6 


«— 32 3 x8 0=6.0 (T = 243 MeV) 




— 32 3 x6 0=6.0 (T = 324 MeV) 




— 32"'x4 0=6.0 (T = 486 MeV) 






•— 64 x16 0=6.0 (T= 121 MeV) 

- - 64"'xl2 0=6.0 (T = 162 MeV) 

• • 64 3 x8 0=6.0 (T = 243 MeV) 
*— 68 3 x8 0=6.0347 (T = 260 MeV) 

52 3 x6 0=5.8876 (T == 265 MeV) 

• ' 72 3 x8 0=6.0684 (T = 275 MeV) 
56 3 x6 0=5.9266 (T = 285 MeV) 

• — ' 76"'x8 0=6.1009 (T == 290 MeV) 

- 60 3 x6 0=5.9640 (T == 305 MeV) 

- ■ 80 3 x8 p=6.1326 (T = 305 MeV) 
— • 64 3 x6 p=6.0 (T = 324 MeV) 
— ■ 64 3 x4 p=6.0 (T = 486 MeV) 



Figure 4: Transverse gluon propagator for ~ (3.3fm) 3 (left) and ~ (6.5fm) 3 (right) spatial lattice volumes. 

The reader should note that the simulation parameters have been carefully chosen, such that 
we only consider two different spatial physical volumes: ~ (3.3fm) 3 and ~ (6.5fm) 3 . This allows 
for a better control of finite size effects. 

In figures [4| and || we show our results. In what concerns the transverse propagator, it decreases 
with the temperature for low momenta. Moreover, this component shows finite volume effects; in 
particular, the large volume data exhibits a turnover in the infrared, not seen at the small volume 
data. The longitudinal component increases for temperatures below T c ~ 270 MeV. Around T Ci 
there is a discontinuity and the propagator decreases for T > T c . The behaviour of the gluon 
propagator as a function of the temperature can also be seen in the 3d plots shown in figure H. 




-•32 x16 0=^.0 (T = 121 MeV) 

- 32 3 xl4 0=6.0 (T = 139 MeV) 

- 32 3 xl2 0=6.0 (T= 162 MeV) 
-* 32*110 0=6.0 (T= 194 MeV) 

- 32 3 x8 0=6.0 (T = 243 MeV) 

- 32 3 x6 0=6.0 (T = 324 MeV) 

- 32 3 x4 $=6.0 (T = 486 MeV) 




— •64 x16 0=6.0 (T=121 McV| 
— ■ 64 3 xl2 0=6.0 (T = 162 MeV) 

— 64 3 x8 0=6.0 (T = 243 MeV ) 
— * 68 3 x8 0=6.0347 (T = 260 MeV) 

52 3 x6 0=5.8876 (T = 265 MeV) 
— • 72 3 x8 0=6.0684 (T = 275 MeV) 

56*6 0=5.9266 (T = 285 MeV) 
< — ' 76 3 x8 0=6.1009 (T = 290 MeV) 

— 60 3 x6 0=5.9640 (T = 305 MeV) 

— 80 3 x6 0=6.1326 (T = 305 MeV) 
— • 64 3 x6 0=6.0 (T = 324 MeV ) 
— ■ 64 3 x4 0=6.0 (T = 486 MeV) 



Figure5: Longitudinal gluon propagator for ~ (3.3fm) 3 (left)and^ (6.5fm) 3 (right) spatial lattice volumes. 





Figure 6: Longitudinal (left) and transverse (right) gluon propagator as a function of momentum and tem- 
perature for a ~ (6.5fm) 3 spatial lattice volume. 
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The reported finite volume effects can be seen in more detail in figure [7] for a temperature of 
324 MeV. 

A complete analysis should include a study of finite lattice spacing effects. For such a goal, we 
worked out two different simulations, at T=305 MeV, with similar physical volumes but different 
lattice spacings. Figure || shows that, for this temperature, finite lattice spacing effects seem to be 
under control, with the exception of the zero momentum for the transverse component. 





Figure 8: Longitudinal (left) and transverse (right) gluon propagator for different lattice spacings (but simi- 
lar physical volume) at T=305 MeV. 

Our preliminary results at finite temperature are rather similar to previous works — see for 
example [Q] and references therein. We are currently working on a complete analysis of our data, 
towards a better understanding of lattice effects on the gluon propagator at finite temperature. 
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